Ecological analysis of world trade 

L.Ermann 1 and D.L.Shcpelyansky 2 

1 Departamento de Fisica, CNEA, Libertador 8250, (C1429BNP) Buenos Aires, Argentina 
2 Laboratoire de Physique Theorique du CNRS, IRSAMC, Universite de Toulouse, UPS, 31062 Toulouse, France 

(Dated: January 17, 2012) 

PACS numbers: 87.23.-n,89.65.Gh, 89.75.Fb, 89.75.Hc 



Ecological systems have a high level of com- 
plexity combined with stability and rich biodi- 
versity. Recently, the analysis of their properties 
and evolution has been pushed forward on a basis 
of concept of mutualistic networks that provides 
a detailed understanding of their features being 
linked to a high nestedness of these networks. 
It was shown that the nestedness architecture of 
mutualistic networks of plants and their pollina- 
tors minimizes competition and increases biodi- 
versity. Here, using the United Nations COM- 
TRADE database for years 1962 - 2009, we show 
that a similar ecological analysis gives a valuable 
description of the world trade. In fact the coun- 
tries and trade products are analogous to plants 
and pollinators, and the whole trade network is 
characterized by a low nestedness temperature 
which is typical for the ecological networks. This 
approach provides new mutualistic features of the 
world trade highlighting new significance of coun- 
tries and trade products for the world trade. 

Ecological systems are characterized by high complex- 
ity and biodiversity [1] linked to nonlinear dynamics and 
chaos emerging in the process of their evolution [2, 3]. 
The interactions between species form a complex network 
whose properties can be analyzed by the modern meth- 
ods of scale- free networks [4-7] . An important feature of 
ecological networks is that they are highly structured, be- 
ing very different from randomly interacting species [7, 8] . 
Recently is has been shown that the mutualistic networks 
between plants and their pollinators [8-11] are character- 
ized by high nestedness [12, 13] which minimizes compe- 
tition and increases biodiversity. It is argued [12] that 
such type of networks appear in various social contexts 
such as garment industry [13] and banking [14, 15]. Here 
we apply a nestedness analysis to the world trade net- 
work using the United Nations COMTRADE database 
[16] for the years 1962 - 2009. Our analysis shows that 
countries and trade products have relations similar to 
those of plants and pollinators and that the world trade 
network is characterized by a high nestedness typical of 
ecosystems [11, 12]. This provides new mutualistic char- 
acteristics for the world trade. 

Results 

The mutualistic World Trade Network (WTN) is con- 




FIG. 1: Normalized import /export WTN matrix elements 
m' ! ' and m' e ' shown on left /right panels for years 1968 (bot- 
tom) and 2008 (top). Here products/countries (p = 1, . . . , N p 
and c — 1, . . . , N c ) are ordered in a decreasing order of prod- 
uct/country total import or export in a given year. 



structed on the basis of the UN COMTRADE database 
[16] from the matrix of trade transactions MZ expressed 
in USD for a given product (commodity) p from coun- 
try c to country d in a given year (from 1962 to 2009). 
For product classification we use 3-digits Standard In- 
ternational Trade Classification (SITC) Rev.l with the 
number of products N p = 182. All these products are 
described in [16] in the commodity code document SITC 
Revl. The number of countries varies between N c = 164 
in 1962 and iV c = 227 in 2009. The import and ex- 
port trade matrices are defined as M$ c = J2c>=i M L> 
and Mpfj = J2^=i c respectively. We use the di- 
mensionless matrix elements my> = Afw /M max and 
77i( e ) = M^ e > /Mmax where for a given year M max = 
max{max [Mpl] , max [M$]}. The distribution of ma- 
trix elements m^\ in the plane of indexes p and 
c, ordered by the total amount of import/export in a 
decreasing order, are shown in Fig. 1 for years 1968 and 
2008 (years 1978, 1988, 1998 are shown in Fig. S-l of Sup- 
porting Information (SI)). These Figs, show that globally 
the distributions of m' 1 ', m!^ 1 remain stable in time es- 
pecially in a view of 100 times growth of the total trade 
volume during the period 1962-2009. The fluctuations of 
m( e ) arc visibly larger compared to mW case since cer- 
tain products, e.g. petroleum, arc exported by only a few 
countries while it is imported by almost all countries. 

To use the methods of ecological analysis we construct 
the mutualistic network matrix for import and ex- 



port Q( e ' whose matrix elements take binary value 1 or 
if corresponding elements m^' and m/ e ' are respec- 
tively larger or smaller than a certain trade threshold 
value /i. The fraction ip of nonzero matrix elements varies 
smoothly in the range 10~ 6 < fj, < 10~ 2 (see Fig. S-2 of 
SI) and the further analysis is not really sensitive to the 
actual /i value inside this broad range. It is important 
to note that in contrast to ecological systems [12] the 
world trade is described by a directed network and hence 
we characterize the system by two mutualistic matrices 
and corresponding to import and export. Using 
the standard nestedness BINMATNEST algorithm [17] 
we determine the nestedness parameter r\ of the WTN 
and the related nestedness temperature T ~ 100(1 — rj). 
The algorithm reorders lines and columns of a mutual- 
istic matrix concentrating nonzero elements as much as 
possible in the top left corner and thus providing infor- 
mation about causal role of immigration and extinction 
in an ecological system. A high level of nestedness and 
ordering can be reached only for systems with low T . It 
is argued that the nested architecture of real mutualistic 
networks increases their biodiversity. 

The nestedness matrices generated by the BIN- 
MATNEST algorithm [17] arc shown in Fig. 2 for ecol- 
ogy networks ARM (N pl = 84, N amm = 101, if = 0.043, 
T = 2.4) and WES (N pl = 207, N amrn = 110, ip = 0.049, 
T = 3.2) from [11, 18]. Using the same algorithm we 
generate the nestedness matrices of WTN using the mu- 
tualistic matrices for import and export for the 
WTN in years 1968 and 2008 using a fixed typical thresh- 
old fi = 10~ 3 (see Fig. 2; the distributions for other \i- 
values have a similar form and are shown in Fig. S-3 of 
SI). As for ecological systems, for the WTN data we also 
obtain rather small nestedness temperature (T w 6/8 
for import/export in 1968 and T « 4/8 in 2008 respec- 
tively). These values are by a factor 9/4 of times smaller 
than the corresponding T values for import / export from 
random generated networks with the corresponding val- 
ues of ip. The detailed data for T in all years are shown 
in Fig. 3 and the comparison with the data for random 
networks are given in Figs. S-4,S-5,S-6 in SI. The data 
of Fig. 3 show that the value of T changes by about 
30-40% with variation of \i by a factor 1000. We think 
that this is relatively small variation of T compared to 
enormous variation of ji that confirms the stability and 
relevance of ecological analysis and nestedness ordering. 
The nestedness temperature T remains rather stable in 
time: in average there is 40% drop of T from 1962 to 
2000 and 20% growth from 2000 to 2009. We attribute 
the growth in last decade to the globalization of trade. 
The small value of nestedness temperature obtained for 
the WTN confirms the validity of the ecological analysis 
of WTN structure: trade products play the role of polli- 
nators which produce exchange between world countries, 
which play the role of plants. Like in ecology the WTN 
evolves to the state with very low nestedness tempera- 
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FIG. 2: Nestedness matrices for the plant-animal mutualis- 
tic networks on top panels, and for the WTN of countries- 
products on middle and bottom panels. Top-left and top- 
right panels represent data of ARR1 and WES networks from 
[11, 18]. The WTN matrices are computed with the threshold 
H — 10 -3 and corresponding ip ^ 0.2 for years 1968 (bottom) 
and 2008 (middle) for import (left panels) and export (right 
panels). Red and blue represent unit and zero elements re- 
spectively; only lines and columns with nonzero elements are 
shown. The order of plants-animals, countries-products is 
given by the nestedness algorithm [17], the perfect nestedness 
is shown by green curves for the corresponding values of ip, 

ture that satisfies the ecological concept of system stabil- 
ity appearing as a result of high network nestendess [12]. 

The nestedness algorithm [17] creates effective ecologi- 
cal ranking (EcoloRanking) of all UN countries. The evo- 
lution of 20 top ranks throughout the years is shown in 
Fig. 4 for import and export. This ranking is quite differ- 
ent from the more commonly applied ranking of countries 
by their total import/export monetary trade volume [19] 
(see corresponding data in Fig. 5) or recently proposed 
democratic ranking of WTN based on the Google matrix 
analysis [20]. Indeed, in 2008 China is at the top rank 
for total export volume but it is only at 5th position in 
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FIG. 3: Nestedness temperature T as a function of years 
for the WTN for ju = 10" 3 (curves), 10" 4 (circles), 10" 6 
(squares); import and export data are shown in red and blue. 



EcoloRank (see Fig. 4, Fig. 5 and Table I in SI). In a 
similar way Japan moves down from 4th to 17th position 
while USA raises up from 3rd to 1st rank. 

The same nestedness algorithm generates not only the 
ranking of countries but also the ranking of trade prood- 
ucts for import and export which is presented in Fig. 6. 
For comparison we also show there the standard ranking 
of products by their trade volume. In Fig. 6 the color of 
symbol marks the 1st SITC digit described in [16] and in 
Table II in SI. 

Discussion 

The origin of such a difference between EcoloRanking 
and trade volume ranking of countries is related to the 
main idea of mutualistic ranking in ecological systems: 
the nestedness ordering stresses the importance of mu- 
tualistic pollinators (products for WTN) which generate 
links and exchange between plants (countries for WTN). 
In this way generic products, which participate in the 
trade between many countries, become of primary impor- 
tance even if their trade volume is not at the top lines 
of import or export. In fact such mutualistic products 
glue the skeleton of the world trade while the nestedness 
concept allows to rank them in order of their importance. 
The time evolution of this EcoloRanking of products of 
WTN is shown in Fig. 6 for import / export in comparison 
with the product ranking by the monetary trade volume 
(since the trade matrix is diagonal in product index the 
ranking of products in the latter case is the same for im- 
port/export). The top and middle panels have dominate 
colors corresponding to machinery (SITC 7; blue) and 
mineral fuels (3; black) with a moderate contribution of 
chemicals (5; yellow) and manufactured articles (8; cyan) 
and a small fraction of goods classified by material (6; 
green). Even if the global structure of product ranking 
by trade volume has certain similarities with import Ecol- 
oRanking there are also important new elements. Indeed, 
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FIG. 4: Top 20 EcoloRank countries as a function of years for 
the WTN import/export on top/bottom panels. The ranking 
is given by the nestedness algorithm [17] for the trade thresh- 
old fj, = 10~ 3 ; each country is represented by its corresponding 
flag. As an example, dashed lines show time evolution of the 
following countries: USA, UK, Japan, China, Spain. 



in 2008 the mutualistic significance of petroleum prod- 
ucts (SITC 332), machindus (machines for special indus- 
tries 718) and medpharm (mcdical-pharmaceutic prod- 
ucts 541) is much higher compared to their volume rank- 
ing, while petroleum crude (331) and office machines 
(714) have smaller mutualistic significance compared to 
their volume ranking. 

The new element of EcoloRanking is that it differ- 
entiates between import and export products while for 
trade volume they are ranked in the same way. In- 
deed, the dominant colors for export (Fig. 6 bottom 
panel) correspond to food (SITC 0; red) with contri- 
bution of black (present in import) and crude mate- 
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FIG. 5: Top 20 countries as a function of years ranked by the 
total monetary trade volume of the WTN in import /export on 
top/bottom panels respectively; each country is represented 
by its corresponding flag. Dashed lines show time evolution 
of the same countries as in Fig. 4. 
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FIG. 6: Top 10 ranks of trade products as a function of years 
for the WTN. Top panel: ranking of products by monetary 
trade volume; middle/bottom panels: ranking is given by the 
nestedness algorithm [17] for import/export with the trade 
threshold fi = 10~ 3 . Each product is shown by its own symbol 
with short name written at years 1968, 2008; symbol color 
marks 1st SITC digit; SITC codes of products and their names 
are given in Table II of SI. 



rials (2; violet); followed by cyan (present in import) 
and more pronounced presence of finnotclass (commodi- 
ties/transactions not classified 9; brown). EcoloRanking 
of export shows a clear decrease tendency of dominance of 
SITC0 and SITC2 with time and increase of importance 
of SITC3,7. It is interesting to note that petroleum prod- 
ucts SITC332 is vary vulnerable in volume ranking due 
to significant variations of petroleum prices but in Ecol- 
oRanking this product keeps the stable top positions in 
all years showing its mutualistic structural importance 
for the world trade. EoloRanking of export shows also 
importance of fish (SITC031), clothing (SITC841) and 
fruits (SITC051) which are placed on higher positions 
compared to their volume ranking. At the same time 



roadvehic (SITC732), which are at top volume ranking, 
have relatively low ranking in export since only a few 
countries dominate the production of road vehicules. 

It is interesting to note that in Fig. 6 petroleum crude 
is at the top of trade volume ranking e.g. in 2008 (top 
panel) but it is absent in import EcoloRanking (middle 
panel) and it is only on 6th position in export EcoloRank- 
ing (bottom panel). A similar feature is visible for years 
1968, 1978. On a first glance this looks surprising but 
in fact for mutualistic EcoloRanking it is important that 
a given product is imported from top EcoloRank coun- 
tries: this is definitely not the case for petroleum crude 
which practically is not produced inside top 10 import 
EcoloRank countries (the only exception is USA, which 
however also does not export much). Due to that reason 
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this product has low mutualistic significance. 

The mutualistic concept of product importance is at 
the origin of significant difference of EcoloRanking of 
countries compared to the usual trade volume ranking 
(see Fig. 4, Fig. 5). Indeed, in the latter case China and 
Japan arc at the dominant positions but their trade is 
concentrated in specific products which mutualistic role 
is relatively low. In contrast USA, Germany and France 
keep top three EcoloRank positions during almost 40 
years clearly demonstrating their mutualistic power and 
importance for the world trade. 

In conclusion, our results show the universal features 
of ecologic ranking of complex networks with promising 
future applications to trade, finance and other areas. 

Acknowledgments: We thank Arlene Adriano and 
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SUPPORTING INFORMATION 

Here we present the Supporting Information (SI) for 
the main part of the paper. 

In Fig. S-l, in a complement to Fig. 1, we show the nor- 
malized WTN matrix for import mM> and export at 
additional year 1978,1988,1998. As in Fig. 1 all prod- 
ucts and countries are ordered in a decreasing order of 
product (p = 1, . . . , N — p) and country (c = 1, . . . , N c ) 
import (left panels) and export (right panels) in a given 
year. These data show that the global distribution re- 
mains stable in time: indeed, the global monetary trade 
volume was increased by a factor 100 from year 1962 to 
2008 (sec e.g. Fig. 5 in [20]) but the shape of the distri- 
bution remained essentially the same. 




FIG. S-l: Same type of WTN matrix data as in Fig. 1 shown 
for years 1978, 1988, 1998 in panels from bottom to top re- 
spectively. 



The dependence of the fraction <p of nonzero elements 
of the mutualistic matrices of import and export 
on the cutoff threshold fi is shown in Fig. S-2. In 
the range of 10~ 6 < fj, < 10 -2 there is a smooth relatively 
weak variation of ip with p. 

In Fig. S-3, in addition to Fig. 2, we show the nest- 
edness matrices of WTN at various values of the cutoff 
threshold (i. The data at various fi values show that in 
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FIG. S-2: The fraction ip of nonzero matrix elements for the 
mutualistic network matrices of import and export 
as a function of the cutoff trade threshold [i for the normalized 
WTN matrices m*- 1 ' and for the year 2008; the red curve 
shows the case of import while the blue curve shows the case 
of export network. 

all cases the nestedness algorithm [17] correctly generates 
a matrix with nestedness structure. 



imports exports 




countries 



FIG. S-3: Same as in Fig. 2: nestedness matrix for the 
WTN data in 2008 shown for the threshold values /i = 
10 -6 , 10 -4 , 10 -2 ) (from top to bottom); the perfectly nest- 
edness is shown by green curves for the corresponding values 
of tp taken from Fig. S-2. 

The variation of the nestedness temperature T with 
time is shown in Fig. 3 and several values of the trade 



threshold \i. These data show that in average the value 
of T for export is higher than for import. We attribute 
this to stronger fluctuations of matrix elements of 
compared to those of that is well visible in Figs. 1,S- 
1. As it is pointed in the main part, we attribute this to 
the fact that e.g. only a few countries export petroleum 
crude while the great majority of countries import this 
product. 
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year 

FIG. S-4: Nestedness temperature T for the model given by 
random generated networks; here T is computed with 500 ran- 
dom realisations of network for each year using N p , N c and ip 
of the corresponding WTN data in this year at fj, = 10 -3 ; im- 
port/export data are shown by red/blue curves respectively. 

In Fig. S-4 we show the nestedness temperature depen- 
dence on time for the case of random generated networks 
which have the same fraction of nonzero matrix elements 
ip as the WTN at the given year and \i = 10~ 3 . These 
data, compared with those of Fig. 3, really demonstrate 
that the real WTN has values of T by a factor 5 (export) 
to 10 (import) smaller comparing to the random net- 
works. This confirm the nestedness structure of WTN 
being similar to the case of ecology networks discussed 
in [12]. It is interesting to note that for random gener- 
ated networks the values of T for import are larger than 
for export while to the WTN we have the opposite re- 
lation. The histogram of distribution of T for random 
generated networks for all years 1962-2009 is shown in 
Fig. S-5. Even minimal values of T remain several times 
larger than the WTN values of T. 

In Fig. S-6 we show the dependence of T on the trade 
threshold // for the WTN data in year 2008. We see 
that there is only about 10-20% of variation of T for the 
range 10 -5 < fj, < 10~ 3 . Even for a much larger range 
10~ 6 < n < 10~ 2 the variation of T remains smooth 
and remains in the bounds of 100%. This confirms the 
stability of nestedness temperature in respect to broad 
range variations of fi. We present the majority of our 
data for fj, = 10~ 3 which is approximately located in the 
flat range of T variation in year 2008. The data of Table I 
for EcoloRanking of countries at two different values of 
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FIG. S-5: Histogram of temperatures for 500 random gener- 
ated networks per year (from 1962 to 2009). Top (bottom) 
panel represents import (export) data; here the parameter 
values of N p , N c and cp are as for the corresponding WTN 
years at /i = 10 -3 . 



H in year 2008 confirm the stability of this ncstedncss 
ordering. At the same time larger values of \x stress the 
importance of countries with a large trade volume, e.g. 
the position of China in export goes up from rank 5 at 
fj, = 10~ 3 to rank 3 at fx = 10~ 2 . 



In Table I we present trade volume ranking and Ecol- 
oRanking of top 20 countries for import/export of WTN 
in year 2008. 



T 

10 
5 




qI I 

io" 6 io" 5 10" 4 II io" 3 10" 2 10"' 



FIG. S-6: Nestedness temperature in the WTN at year 2008 
as a function of threshold n; imports/exports networks are 
shown by red/blue curves respectively. 



In Table II we give the notations and symbols for Fig. 6 
with corresponding SITC Revl codes and names. The 
list of all SITC Revl codes is available at [16] (sec file 

http : //unstats ,un. org/unsd/tradekb/Attachmentl93 . aspx). 

The colors of symbols in Fig. 4 mark the first digit of 
SITC Revl code: - red (Food and live animals); 1 

- does not appear in Fig. 4 (Beverages and tobacco); 
2 - violet (Crude materials, inedible, except fuels); 3 - 
black (Mineral fuels, lubricants and related materials); 4 

- does not appear in Fig. 4 (Animal and vegetable oils 
and fats); 5 - yellow (Chemicals); 6 - green (Manufact 
goods classified chiefly by material) ; 7 - blue (Machinery 
and transport equipment); 8 - cyan (Miscellaneous 
manufactured articles); 9 - brown (Commod. and 
transacts. Not class. Accord. To kind). 
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TABLE I: Top 20 ranks of countries for import and export with ranking by the monetary trade volume and by the nestedness 
algorithm at two threshold values /x (year 2008). 
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Netherlands 


tnma 


UK 


Spain 


Japan 


1U 


Canada 


bnina 


Spain 


Belgium 


Canada 


Spain 


i i 
11 


Spain 


Canada 


Canada 


Canada 


India 


Canada 


1 o 
Iz 


Rep. of Korea 


Mexico 


Russian Fed. 


Rep. of Korea 


Poland 


Switzerland 


1 Q 

lo 


Russian Fed. 


Rep. of Korea 


Rep. of Korea 


Mexico 


Sweden 


India 


1 A 
14 


Mexico 


Russian Fed. 


Switzerland 


Saudi Arab 


Austria 


Rep. of Korea 


15 


Singapore 


Poland 


Austria 


Singapore 


Brazil 


Poland 


16 


India 


Austria 


Poland 


Spain 


Australia 


Turkey 


17 


Poland 


Switzerland 


Sweden 


Malaysia 


Japan 


Czech Rep. 


18 


Switzerland 


Turkey 


Mexico 


Brazil 


Russian Fed. 


Austria 


19 


Turkey 


U. Arab Emir. 


India 


India 


Denmark 


Thailand 


20 


Brazil 


Denmark 


Singapore 


Switzerland 


Thailand 


Denmark 
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TABLE II: Product names for SITC Revl 3-digit code used in Fig. 4 



Symbol Code Abbreviation Name 





001 


animals 


Live animals 




031 


fish 


Fish, fresh & simply preserved 




051 


fruits 


Fruit, fresh, and nuts excl. Oil nuts 




054 


vegetables 


Vegetables, roots & tubers, fresh or dried 




061 


sugarhon 


Sugar and honey 




071 


coffee 


Coffee 




081 


feedanim 


Feed, stuff for animals excl. unmilled cereals 




221 


oilseeds 


Oil seeds, oil nuts and oil kernels 




263 


cotton 


Cotton 




283 


ores 


Ores & concentrates of non ferrous base metals 




331 


petrolcrude 


Petroleum, crude and partly refined 




332 


petrolprod 


Petroleum products 




341 


gas 


Gas, natural and manufactured 




512 


orgchem 


Organic chemicals 




541 


medpharm 


Medicinal & pharmaceutical products 




581 


plasticmat 


Plastic materials,regenerd. cellulose & resins 




599 


chemmat 


Chemical materials and products, nes 


• 


652 


cottwoven 


Cotton fabrics, woven ex. narrow or spec. fabrics 


■ 


653 


ncottwov 


Text fabrics woven ex narrow, spec, not cotton 


♦ 


667 


pearlsprec 


Pearls and precious and semi precious stones 


A 


674 


iron 


Universals, plates and sheets of iron or steel 


-< 


682 


copper 


Copper 


• 


711 


nelecmach 


Power generating machinery, other than electric 


■ 


714 


offmach 


Office machines 


♦ 


718 


machindus 


Machines for special industries 


▲ 


719 


mapplpart 


Machinery and appliances non electrical parts 


-< 


722 


elecmach 


Electric power machinery and switchgear 


T 


724 


telecomm 


Telecommunications apparatus 


► 


729 


oelecmach 


Other electrical machinery and apparatus 


+ 


732 


roadvehicles 


Road motor vehicles 


X 


735 


ships 


Ships and boats 


• 


841 


clothing 


Clothing except fur clothing 


• 


931 


finnotclass 


Special transactions not classd. accord. to kind 



